<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<!-- 
	Copyright (C) 2007, 2008, 2009, 2010, 2011. PARP Research Group.
	<http://perception.inf.um.es>
	University of Murcia, Spain.

	This file is part of the QVision library.

	QVision is free software: you can redistribute it and/or modify
	it under the terms of the GNU Lesser General Public License as
	published by the Free Software Foundation, version 3 of the License.

	QVision is distributed in the hope that it will be useful,
	but WITHOUT ANY WARRANTY; without even the implied warranty of
	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
	GNU Lesser General Public License for more details.

	You should have received a copy of the GNU Lesser General Public
	License along with QVision. If not, see <http://www.gnu.org/licenses/>.
-->

<html><head><meta http-equiv="content-Type" content="text/html;charset=UTF-8">
<title>QVision: Qt&#39;s Image, Video and Computer Vision Library</title>
<meta name="title" content="QVision" />
<meta name="dc.title" content="QVision" />
<meta name="url" content="http://perception.inf.um.es/QVision" />
<meta name="author" content="PARP Research Group - http://perception.inf.um.es" />
<meta name="revisit-after" content="30 DAYS"/>
<meta name="robots" content="index,follow"/>
<meta name="classification" content="*">
<meta name="rating" content="Safe For Kids">
<meta name="distribution" content="GLOBAL"/>
<meta name="description" content="Qt's Image, Video and Computer Vision Library"/>
<meta name="page-topic" content="Computer Vision research and prototype programming"/>
<meta name="geo.country" content="ES" />

<!--
Keywords:
By license:		open source, gnu, lgpl, gpl, free
By theme:		computer vision, image processing, robotics, programming, source, development
By usage:		library, toolkit, framework, prototype, application
By programming specs:	object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping
Interoperability with:	Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack
Functionallity:		image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface
Main data-types:	matrix, vector, tensor, quaternion, image, polyline
Video sources:		webcam, camera, stream
Devices:		embedded, desktop computer, laptop, mini-laptop
Authors:		PARP research group. University of Murcia, Spain.
-->

<meta name="keywords" content="augmented reality, sfm, structure from motion, open source, gnu, lgpl, gpl, free, computer vision, image processing, robotics, programming, source, development, library, toolkit, framework, prototype, application, object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping, Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack, image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface, matrix, vector, tensor, quaternion, image, polyline, webcam, camera, stream, embedded, desktop computer, laptop, mini-laptop, University of Murcia, Spain, PARP research group, vision por computador"/>
<meta http-equiv="keywords" content="augmented reality, sfm, structure from motion, open source, gnu, lgpl, gpl, free, computer vision, image processing, robotics, programming, source, development, library, toolkit, framework, prototype, application, object oriented, c++, block programming, reusability, gui, graphical, parallel computing, high performance, GPU, prototyping, Qt, GSL, GNU Scientific library, OpenCV, CGAL, QWT, CUDA, mplayer, IPP, Intel Image Performance Primitives, blas, lapack, image features, matrix algebra, projective geometry, mser, function minimization, function optimization, canny operator, harris operator, corner detection, performance evaluation, cpu usage, graphical interface, matrix, vector, tensor, quaternion, image, polyline, webcam, camera, stream, embedded, desktop computer, laptop, mini-laptop, University of Murcia, Spain, PARP research group, vision por computador"/>
<meta http-equiv="pragma" content="no-cache"/>
<meta http-equiv="title" content="QVision"/>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="tabs.css" rel="stylesheet" type="text/css" />
<link rel="shortcut icon" href="favicon.ico" />
</head><body>

<table width="100%"><tr>
	<td><a href="http://perception.inf.um.es/"><img src="parp.png" border="0" /> <big>PARP Research Group</big></a></td>
	<td align="right"><a href="http://www.um.es/"><big>Universidad de Murcia</big> <img src="um.png" border="0" /></a></td>
</tr></table>

<hr /><br />

<table width="95%" align="center"><tr><td>

<!-- Generated by Doxygen 1.6.3 -->
<script type="text/javascript"><!--
var searchBox = new SearchBox("searchBox", "search",false,'Search');
--></script>
<div class="navigation" id="top">
  <div class="tabs">
    <ul>
      <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
      <li><a href="modules.html"><span>Modules</span></a></li>
      <li><a href="annotated.html"><span>Classes</span></a></li>
      <li class="current"><a href="files.html"><span>Files</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
          <form id="FSearchBox" action="search.php" method="get">
            <img id="MSearchSelect" src="search/search.png" alt=""/>
            <input type="text" id="MSearchField" name="query" value="Search" size="20" accesskey="S" 
                   onfocus="searchBox.OnSearchFieldFocus(true)" 
                   onblur="searchBox.OnSearchFieldFocus(false)"/>
          </form>
        </div>
      </li>
    </ul>
  </div>
  <div class="tabs">
    <ul>
      <li><a href="files.html"><span>File&nbsp;List</span></a></li>
      <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
    </ul>
  </div>
<h1>examples/OpenCV/mser_sample/mser_sample.cpp</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/* This sample code was originally provided by Liu Liu</span>
<a name="l00002"></a>00002 <span class="comment"> * Copyright� 2009, Liu Liu All rights reserved.</span>
<a name="l00003"></a>00003 <span class="comment"> */</span>
<a name="l00004"></a>00004 
<a name="l00005"></a>00005 <span class="comment">/*</span>
<a name="l00006"></a>00006 <span class="comment">This source is one of the sample codes from the OpenCV 2.3 library</span>
<a name="l00007"></a>00007 <span class="comment">*/</span>
<a name="l00008"></a>00008 
<a name="l00009"></a>00009 <span class="preprocessor">#include &quot;opencv2/highgui/highgui.hpp&quot;</span>
<a name="l00010"></a>00010 <span class="preprocessor">#include &quot;opencv2/features2d/features2d.hpp&quot;</span>
<a name="l00011"></a>00011 <span class="preprocessor">#include &quot;opencv2/imgproc/imgproc_c.h&quot;</span>
<a name="l00012"></a>00012 
<a name="l00013"></a>00013 
<a name="l00014"></a>00014 <span class="keywordtype">void</span> help()
<a name="l00015"></a>00015 {
<a name="l00016"></a>00016         printf(<span class="stringliteral">&quot;\nThis program demonstrates the Maximal Extremal Region interest point detector.\n&quot;</span>
<a name="l00017"></a>00017                         <span class="stringliteral">&quot;It finds the most stable (in size) dark and white regions as a threshold is increased.\n&quot;</span>
<a name="l00018"></a>00018                         <span class="stringliteral">&quot;\nCall:\n&quot;</span>
<a name="l00019"></a>00019                         <span class="stringliteral">&quot;./mser_sample &lt;path_and_image_filename, Default is &#39;puzzle.png&#39;&gt;\n\n&quot;</span>);
<a name="l00020"></a>00020 }
<a name="l00021"></a>00021 
<a name="l00022"></a>00022 <span class="keyword">static</span> CvScalar colors[] =
<a name="l00023"></a>00023 {
<a name="l00024"></a>00024     {{0,0,255}},
<a name="l00025"></a>00025     {{0,128,255}},
<a name="l00026"></a>00026     {{0,255,255}},
<a name="l00027"></a>00027     {{0,255,0}},
<a name="l00028"></a>00028     {{255,128,0}},
<a name="l00029"></a>00029     {{255,255,0}},
<a name="l00030"></a>00030     {{255,0,0}},
<a name="l00031"></a>00031     {{255,0,255}},
<a name="l00032"></a>00032     {{255,255,255}},
<a name="l00033"></a>00033     {{196,255,255}},
<a name="l00034"></a>00034     {{255,255,196}}
<a name="l00035"></a>00035 };
<a name="l00036"></a>00036 
<a name="l00037"></a>00037 <span class="keyword">static</span> uchar bcolors[][3] =
<a name="l00038"></a>00038 {
<a name="l00039"></a>00039     {0,0,255},
<a name="l00040"></a>00040     {0,128,255},
<a name="l00041"></a>00041     {0,255,255},
<a name="l00042"></a>00042     {0,255,0},
<a name="l00043"></a>00043     {255,128,0},
<a name="l00044"></a>00044     {255,255,0},
<a name="l00045"></a>00045     {255,0,0},
<a name="l00046"></a>00046     {255,0,255},
<a name="l00047"></a>00047     {255,255,255}
<a name="l00048"></a>00048 };
<a name="l00049"></a>00049 
<a name="l00050"></a>00050 
<a name="l00051"></a>00051 <span class="keywordtype">int</span> main( <span class="keywordtype">int</span> argc, <span class="keywordtype">char</span>** argv )
<a name="l00052"></a>00052 {
<a name="l00053"></a>00053         <span class="keywordtype">char</span> path[1024];
<a name="l00054"></a>00054         IplImage* img;
<a name="l00055"></a>00055         help();
<a name="l00056"></a>00056         <span class="keywordflow">if</span> (argc!=2)
<a name="l00057"></a>00057         {
<a name="l00058"></a>00058                 strcpy(path,<span class="stringliteral">&quot;puzzle.png&quot;</span>);
<a name="l00059"></a>00059                 img = cvLoadImage( path, CV_LOAD_IMAGE_GRAYSCALE );
<a name="l00060"></a>00060                 <span class="keywordflow">if</span> (!img)
<a name="l00061"></a>00061                 {
<a name="l00062"></a>00062                         printf(<span class="stringliteral">&quot;\nUsage: mser_sample &lt;path_to_image&gt;\n&quot;</span>);
<a name="l00063"></a>00063                         <span class="keywordflow">return</span> 0;
<a name="l00064"></a>00064                 }
<a name="l00065"></a>00065         }
<a name="l00066"></a>00066         <span class="keywordflow">else</span>
<a name="l00067"></a>00067         {
<a name="l00068"></a>00068                 strcpy(path,argv[1]);
<a name="l00069"></a>00069                 img = cvLoadImage( path, CV_LOAD_IMAGE_GRAYSCALE );
<a name="l00070"></a>00070         }
<a name="l00071"></a>00071 
<a name="l00072"></a>00072         <span class="keywordflow">if</span> (!img)
<a name="l00073"></a>00073         {
<a name="l00074"></a>00074                 printf(<span class="stringliteral">&quot;Unable to load image %s\n&quot;</span>,path);
<a name="l00075"></a>00075                 <span class="keywordflow">return</span> 0;
<a name="l00076"></a>00076         }
<a name="l00077"></a>00077         IplImage* rsp = cvLoadImage( path, CV_LOAD_IMAGE_COLOR );
<a name="l00078"></a>00078         IplImage* ellipses = cvCloneImage(rsp);
<a name="l00079"></a>00079         cvCvtColor(img,ellipses,CV_GRAY2BGR);
<a name="l00080"></a>00080         CvSeq* contours;
<a name="l00081"></a>00081         CvMemStorage* storage= cvCreateMemStorage();
<a name="l00082"></a>00082         IplImage* hsv = cvCreateImage( cvGetSize( rsp ), IPL_DEPTH_8U, 3 );
<a name="l00083"></a>00083         cvCvtColor( rsp, hsv, CV_BGR2YCrCb );
<a name="l00084"></a>00084         CvMSERParams params = cvMSERParams();<span class="comment">//cvMSERParams( 5, 60, cvRound(.2*img-&gt;width*img-&gt;height), .25, .2 );</span>
<a name="l00085"></a>00085 
<a name="l00086"></a>00086         <span class="keywordtype">double</span> t = (double)cvGetTickCount();
<a name="l00087"></a>00087         cvExtractMSER( hsv, NULL, &amp;contours, storage, params );
<a name="l00088"></a>00088         t = cvGetTickCount() - t;
<a name="l00089"></a>00089         printf( <span class="stringliteral">&quot;MSER extracted %d contours in %g ms.\n&quot;</span>, contours-&gt;total, t/((<span class="keywordtype">double</span>)cvGetTickFrequency()*1000.) );
<a name="l00090"></a>00090         uchar* rsptr = (uchar*)rsp-&gt;imageData;
<a name="l00091"></a>00091         <span class="comment">// draw mser with different color</span>
<a name="l00092"></a>00092         for ( <span class="keywordtype">int</span> i = contours-&gt;total-1; i &gt;= 0; i-- )
<a name="l00093"></a>00093         {
<a name="l00094"></a>00094                 CvSeq* r = *(CvSeq**)cvGetSeqElem( contours, i );
<a name="l00095"></a>00095                 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> j = 0; j &lt; r-&gt;total; j++ )
<a name="l00096"></a>00096                 {
<a name="l00097"></a>00097                         CvPoint* pt = CV_GET_SEQ_ELEM( CvPoint, r, j );
<a name="l00098"></a>00098                         rsptr[pt-&gt;x*3+pt-&gt;y*rsp-&gt;widthStep] = bcolors[i%9][2];
<a name="l00099"></a>00099                         rsptr[pt-&gt;x*3+1+pt-&gt;y*rsp-&gt;widthStep] = bcolors[i%9][1];
<a name="l00100"></a>00100                         rsptr[pt-&gt;x*3+2+pt-&gt;y*rsp-&gt;widthStep] = bcolors[i%9][0];
<a name="l00101"></a>00101                 }
<a name="l00102"></a>00102         }
<a name="l00103"></a>00103         <span class="comment">// find ellipse ( it seems cvfitellipse2 have error or sth?</span>
<a name="l00104"></a>00104         <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i = 0; i &lt; contours-&gt;total; i++ )
<a name="l00105"></a>00105         {
<a name="l00106"></a>00106                 CvContour* r = *(CvContour**)cvGetSeqElem( contours, i );
<a name="l00107"></a>00107                 CvBox2D box = cvFitEllipse2( r );
<a name="l00108"></a>00108                 box.angle=(float)CV_PI/2-box.angle;
<a name="l00109"></a>00109 
<a name="l00110"></a>00110                 if ( r-&gt;color &gt; 0 )
<a name="l00111"></a>00111                         cvEllipseBox( ellipses, box, colors[9], 2 );
<a name="l00112"></a>00112                 <span class="keywordflow">else</span>
<a name="l00113"></a>00113                         cvEllipseBox( ellipses, box, colors[2], 2 );
<a name="l00114"></a>00114 
<a name="l00115"></a>00115         }
<a name="l00116"></a>00116 
<a name="l00117"></a>00117         cvSaveImage( <span class="stringliteral">&quot;rsp.png&quot;</span>, rsp );
<a name="l00118"></a>00118 
<a name="l00119"></a>00119         cvNamedWindow( <span class="stringliteral">&quot;original&quot;</span>, 0 );
<a name="l00120"></a>00120         cvShowImage( <span class="stringliteral">&quot;original&quot;</span>, img );
<a name="l00121"></a>00121 
<a name="l00122"></a>00122         cvNamedWindow( <span class="stringliteral">&quot;response&quot;</span>, 0 );
<a name="l00123"></a>00123         cvShowImage( <span class="stringliteral">&quot;response&quot;</span>, rsp );
<a name="l00124"></a>00124 
<a name="l00125"></a>00125         cvNamedWindow( <span class="stringliteral">&quot;ellipses&quot;</span>, 0 );
<a name="l00126"></a>00126         cvShowImage( <span class="stringliteral">&quot;ellipses&quot;</span>, ellipses );
<a name="l00127"></a>00127 
<a name="l00128"></a>00128         cvWaitKey(0);
<a name="l00129"></a>00129 
<a name="l00130"></a>00130         cvDestroyWindow( <span class="stringliteral">&quot;original&quot;</span> );
<a name="l00131"></a>00131         cvDestroyWindow( <span class="stringliteral">&quot;response&quot;</span> );
<a name="l00132"></a>00132         cvDestroyWindow( <span class="stringliteral">&quot;ellipses&quot;</span> );
<a name="l00133"></a>00133         cvReleaseImage(&amp;rsp);
<a name="l00134"></a>00134         cvReleaseImage(&amp;img);
<a name="l00135"></a>00135         cvReleaseImage(&amp;ellipses);
<a name="l00136"></a>00136 
<a name="l00137"></a>00137 }
</pre></div></div>
</td></tr></table>

<br /><hr><br />
<center><a href="http://perception.inf.um.es/QVision">QVision framework</a>.
<a href="http://perception.inf.um.es">PARP research group</a>.
Copyright &copy; 2007, 2008, 2009, 2010, 2011.</center>
<br />
</body>
</html>
